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, Abstract. I investigate through simulations the redshift dependence of 

several lensing measures for two cosmological models, a flat universe with 
a cosmological constant (ACDM), and an open universe (OCDM). I argue 
that quintessence models can be seen as an intermediate model between 
these two models under the weak gravitational lensing perspective. I cal- 
OO I culate for the convergence field the angular power spectrum, the variance 

and higher order moments, the probability distribution function, and the 
Minkowski functionals. I find that in the redshift range examined (from 
^ ■ z = 0.7 to z = 3 various statistics show an increase of the non-Gaussianity 

CN i of the lensing maps for closer sources. Lensing surveys of higher redshift 

■ are expected to be less non-Gaussian than shallow ones because of the 

i non-linear evolution of the density field, and the central limit theorem. 

o 
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A lensed image contains information about the positions of the lens and source 
^ ■ object, and about the lens itself (e.g. its refractive index). In the cosmological 

^ ', context the mass inhomogeneities in the universe act as gravitational lenses to 

source objects far away. The dominant cosmological picture today describes a 
universe that is expanding and has growing density perturbations. The exact 
characterization and quantification of these phenomena are essential pursuits of 
observational Cosmology, and their description by a theoretical model is a major 
' objective for cosmologists. 

Weak gravitational lensing studies can contribute to this enterprise because 
lensing maps depend on the evolution of the large-scale structure and on the 
geometrical properties of the universe (the cosmological expansion). 

In this section I briefly review some of the formalism of cosmic expansion 
and structure formation, and point how these two aspects of cosmic evolution 
are captured by weak gravitational lensing maps. 

1.1. Cosmological Geometry 

In a Friedmann-Robertson- Walker model, assuming isotropy and homogeneity, 
the geometry of the universe can be described by the metric 



dT^ = dt' - a\t) dr^ + ff^{r)d(t)' , (1) 

where t, r and (j) are the comoving coordinates, and fx is the curvature-dependent 
comoving angular distance (= r for a flat universe, = sinhr for open, and = sinr 
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for closed). The expansion factor a{t) is determined by the equation 

a' - = -k , (2) 

in a universe of curvature k (= 0, —1, 1 in the case that it is flat, open, or closed, 

respectively), and of density p. In terms of the density parameters Q\, Cl^n and 
Qr, for a cosmological constant, matter, and radiation, respectively, 

P = ^ (j^A + n^a-^ + nra-^) . (3) 

It is also useful to define the Hubble parameter 

H ^ (4) 

a 

1 /9 

H{a) = Ho[nA{l-a-'^) + nm{a-^ -a-'^)+a-^] , (5) 

and to remember the relation a = l/(z + 1). 

In a quintessence model the cosmological constant is substituted by a dark 
energy component, which has negative pressure, i.e. a negative wq in its equa- 
tion of state pq = wqpq. In the case in which wq is constant the Hubble 
expansion rate is them 

H{a) = Ho [0A(a-3("'«+i) - a'^) + ^^(a-^ - a'^) + a'^] ^''^ , (6) 

which reduces to the cosmological constant case [equation (5)] when wq = — 1. 
1.2. Density Perturbations 

The current standard explanation for the existence of a large-scale structure 
in the matter distribution in the universe is that initial small fluctuations in a 
remote past evolved (grew) due to gravitational attraction. That is the gravi- 
tational instability scenario. The origin and characteristics of the initial inho- 
mogeneities is subject of great interest for Cosmology and fundamental Physics, 
but will not be discussed here. 

The matter distribution can be described by the density contrast 5 = p/p—1, 
where p is the average density, and its evolution in the linear regime is given by 

'5 + 2-5 = At:Gp5 (7) 
a 

The growing mode solution for this equation can be factored into spatial and 
temporal functions 

J(x,i) = D(t)5(x,t„), (8) 

where the growth function is 

Dia) oc- r -da (9) 
a Jo a 
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which can be approximated, apart from a normahzation factor, by (Carroll, 
Press & Turner 1992; Lahav et al. 1991) 



2 V 2 



Hr 



^m{a) = — ( 

H 



r2A(a) = J^A 



,(10) 

(11) 
(12) 



The growth function must be normalized such that D{a = 1) = 1. For an 
Einstein-de Sitter Universe there is no suppression of gravitational clustering, 
and the growth function reduces to the scale factor, D[a) = a. 

For a quintessence model (Ma et al. 1999) ^^(a) = VLm/[Q.m + (1 — 
^^m)^"^"'''], and the growth factor can be approximated by (apart from a nor- 
malization factor) 



DQ{a) 



~ {-WQfa-Q.rn{a) 

X |j^^/7(a) - 1 + J^„(a) + 



(13) 



1 + 



1 + 



-(0.255+0.305t(;Q- 



0.0027 



WQ 



) [1 - (a)] - (0.366 + 0.266w;q • 



1 - ^m{a) 
70 

0.07. 



WQ 



InOm(a) • 
(14) 

I use the expressions (14) and (10), properly normalized, to plot the curves 
in Figure 1 that show the linear growth factor for several cosmological models. 
Note that the effect of a quintessence description for the dark energy component 
is equivalent to an intermediate suppression of growth between a flat universe 
with cosmological constant and an open universe. 

The power spectrum of the matter density perturbation is a powerful sta- 
tistical descriptor of the large-scale structure of the universe. From linear theory 
(which is valid at the largest scales) the power can be written as 



P{k,a) = Ak'^T'^{k)D'^ 



(15) 



where ^ is a normalization, n is the initial power spectrum index, and T{k) is the 
transfer function, which modifies a initial power-law spectrum of density pertur- 
bations, building in information about the energy components of the universe. 
However, linear theory becomes a bad approximation when 5(x, i) 3> 1, which 
happens at small scales at late stages of the evolution of the density fields. The 
non-linear evolution causes the power spectrum at large wavcnumber k to be 
underestimated by linear theory. The corrected calculation of the density field 
evolution in the non-linear regime requires the use of N-body simulations. From 
that one can obtain approximations for the non-linear power spectrum using 
prescriptions such the ones developed by Hamilton et al. (1991), and Scranton 
k. Dodelson (2000). See figure 2 for plots of the non-linear power spectrum of 
quintessence models with various equation-of-state coefficients, and of an open 
universe, with no dark energy. 
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Figure 1. Linear growth of factor. SCDM has = Ij all other 
models have = 0.3, wq is the equation-of-state coefficient for a 
quintessence field model with Qq = 0.7. 




lo" 10"^ 10"' 10° io' 

k(h"'Mpcr' 



Figure 2. Power spectrum of matter density fluctuations. All mod- 
els (quintessence models and OCDM) were normalized to match the 
density anisotropies observed in the cosmic background radiation (Bun 
& White 1997 and Ma et al. 1999). The case wq = — 1 is equivalent 
to a ACDM model. 
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Figure 3. Geometrical weighting of density perturbations in weak 
lensing. Lines labeled by wq represent quintessence models with such 
equation-of-state coefficient. 



1.3. Weak Lensing Dependence on Cosmological Evolution 

The distortion matrix from which the shear and the convergence fields are defined 
is explicitly dependent on the cosmic evolution from the redshift of the source 
Zs until today {z = 0) 

Aij{e, Zs) = -2 r g{z, zs)didj<^{0, z)dz . (16) 

JO 

This expression integrates the evolution in the geometry of the universe, which 
is encoded in the factor g, and the evolution of the gravitational potential ^. 

For the convergence this dependence is represented thought the comoving 
distance ^ ^ 

^(") = / ,2U( f\/U ' 

J a a''^H[a')/Ho 

in the expression 

K{e,Zs) = ^—Um d[9,w )dw , (18) 

2 Jo a(w') 

where 

/ / X fK{w')fK{w -w') 

9{w , w) = -— . (19) 

jKiw) 

Figure 3 shows the geometrical weighting of the density fiuctuations for a 
source at a redshift z given by the factor g/a at equation (18). The maximum 
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Figure 4. 2D power spectrum of lens-planes. 



contribution for lensing comes from structures with redshifts close to half of the 
source redshift. Quintessence models have an intermediate weighting between a 
ACDM and a OCDM model, and the farther the source the more significant is 
the difference between the models. 



2. Lens Evolution 

In this section I use a set of statistics to characterize the evolution of the mat- 
ter density field, or more properly, the evolution of lens-planes. I analyze the 
projected density contrast field of boxes of 100 h^^Mpc of side that evolve in 
time. The simulations use 128^ dark-matter particles (no baryonic matter) in 
two models, ACDM and OCDM. 

Figure 4 shows the evolution of the two-dimensional power spectrum of the 
lenses, and illustrates the growth of structure in the matter density field. It 
is possible to note that small-scale structures (represented by large k) have a 
stronger increase than large-scale structures (small k), as is expected from non- 
linear evolution. The difference between ACDM and OCDM is more notable at 
large redshifts, because the simulations were normalized to the cluster abundance 
today. 

The growth of structure in time can be also appreciated by the increase of 
the field variance, skewness, and kurtosis - see Figure 5. These two last quantities 
quantify the growth of non-Gaussianities in the density field, which in this case 
is originally Gaussian, due to non-linear evolution. The same information can 
be obtained from the PDF (figure 6). 
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Figure 5. The rms, skewness, and kurtosis of the projected density 
contrast field as redshift (left) and scale factor (right) functions for 
ACDM (solid lines) and OCDM (dashed lines). Projected boxes have 
100 Ih-^Mpc, and were smoothed at 1 Ih ^Mpc scale. 




Figure 6. Probability distribution function for the density contrast 
of lens planes with lh~^Mpc smoothing. 
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The last aspect of the density field to be analyzed is its morphology. I use 
the second (contour length) and third (genus) Minkowski functionals parame- 
terized by the first functional (fractional area) to characterize the morphology 
of the lens-planes. Figure 7 shows that up to a redshift z=5 there is little change 
in these two functionals. The third functional is not symmetric in relation to 
the midpoint along the axes of the plot, which again shows the non-Gaussianity 
of the density field. Both models are almost indistinguishable from one another 
under this statistics. The meaning of these observations is that in both models 
the non-trivial morphology, i.e. the morphology distinct from that of a Gaussian 
field, emerges at earlier times than z=5, and that in the late universe the mor- 
phology of the LSS as described by the Minkowski functionals evolves little if the 
functionals are parameterized by the fractional area, making them independent 
of the PDF. 



3. Redshift Dependence of Lensing Measures 

I generated simulated convergence maps using the multiple-plan lens method 
(see Guimaraes, 2002 for more details) for ACDM and OCDM. I used the Hydra 
N-body code (Couchman, Thomas & Pearce, 1995) to evolve 128^ particles from 
a redshift z = 50 to 2 = 0, simulating the large-scale mass distribution in 
the two models. For the ACDM = 0.3, JIa = 0.7) simulation I adopt a 
cluster normalization equivalent to iTg = 0.9, and a simulation box of 141.3 
h~^Mpc. For the OCDM (fi^ = 0.3, JIa = 0.7) simulation I adopt a cluster 
normalization equivalent to erg = 0.84, and a simulation box of 128 h~^Mpc. A 
Hubble parameter h = 0.7 was used in both simulations. 

Another examined model mimics an universe with the geometry of the 
ACDM model but no structure evolution. In this toy-model the structre at 
any redshift has the same statistical properties that today. 

Figure 8 shows a realization of the convergence map for both models with 
a source plane at three redshifts. The redshift dependence is very clear, but 
the differentiation between ACDM and OCDM is less evident. A quantitative 
characterization of the models requires the use of several statistics applied to 
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Figure 8. Convergence maps for ACDM (left) and OCDM (right) for 
source plane at redshifts 0.7, 1.5 and 3 (top to bottom). All maps are 
1.7x1.7 degrees^, are smoothed with a 1 arcmin top-hat radius, and 
use the same color scale (from deep blue for the minimum value to red 
to the maximum value) 



10 Antonio C. C. Guimaraes 




Figure 9. Angular power spectrum of the convergence field for 
ACDM (left) and OCDM (right) for source plane at redshifts 3, 2, 
1.5, 1, 0.7 (top to bottom curves). No smoothing was used. 



many realizations of the lensing maps. I calculate as a function of source redshift 
the following statistics for the convergence: the angular power spectrum, the rms 
(root mean square), the skewness {S = {K^)/a^), the kurtosis {K = {{K^)/cr'^) — 
3), the probability distribution function, and the Minkowski functionals. The 
results shown in the figures are the result of 20 realizations of the convergence 
map for each model, and for each redshift. The variance of the measures are 
shown only for ACDM and serves as a reference for maps of same size of other 
models. The sizes of the used maps are (side of square), 4.6° for z = 0.7, 3.5° 
for z = 1, 2.5° for z = 1.5, 2.2° for z = 2, and 1.8° for z = 3. 

Figure 9 shows the results for the angular power spectrum of the convergence 
field. There is an overall increase of the spectrum amplitude as the redshift of 
the source increases. This is also reflected in the rms plot (Figure 10), since 
the field variance is proportional to the convergence angular power spectrum. 
The same figure shows a decreasing skewness and kurtosis for larger redshifts. 
These two statistics can be considered measures of the non-Gaussianity of the 
field since for a Gaussian field they are zero. 

A power law fit to the rms results yield the following redshift dependence 
for the shear variance 



(20) 
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Figure 10. The rms, skewness, and kurtosis of the convergence field 
as redshift functions for ACDM (soUd Unes) and OCDM (dashed hnes). 
Maps were smoothed at 1 arcmin scale. 



where a = 1.34 ± 0.08 for OCDM, and a = 1.5 ± 0.1 for ACDM. This later 
result coincides with the value obtained by Jain & Seljak (1997), a = 1.52, but 
is lower than the value obtained by Barber (2002), a = 2.07 it 0.04. 

Departures from Gaussianity can also be observed in the probability distri- 
bution function of the lensing maps - figure 11. I plot the PDF for the Gaus- 
sianized ACDM maps to show more clearly the non-Gaussianity of the maps. 
The convergence fields were Gaussianized by expanding them in Fourier modes, 
and reconstructing the fields with the same mode amplitudes, but random mode 
phases. 

Figure 12 shows the results for the second and third Minkowski function- 
als (contour length and topological genus, respectively), parameterized by the 
first functional (fractional area). I choose not to present the first Minkowski 
functional because it does not carry any extra information in relation to the 
PDF. The parameterization by the fractional area A instead of the map thresh- 
old V = k/cTk is aimed to further decouple the information contained in the 
Minkowski functionals from the PDF. The results for the Gaussianized fields 
compared with the original field results show that the third functional G{A) 
puts in evidence the non-Gaussianity of the maps much more clearly than the 
second functional C{A). Note that the third functional does not show a clear 
redshift dependence under the conditions examined. 
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Figure 11. Probability distribution function of the convergence field 
smoothed at 1 arcmin scale. The dashed lines are the PDF of the 
Gaussianized fields. 




Figure 12. Minkowski functionals parameterized by the fractional 
area A of the convergence field smoothed at 1 arcmin scale. The thin 
lines are for the Gaussianized maps. 
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4. Conclusion 

In this paper I have made predictions for the redshift evolution of statistical mea- 
sures of weak lensing maps. These predictions show a relatively small difference 
between OCDM and ACDM models with adiabatic Gaussian initial conditions, 
and will thus also hold in quintessence models in similar conditions. Given the 
robustness of the results in redshift dependence within the class of models tested, 
I propose the study of the redshift dependence of weak gravitational lensing as a 
test for the stability of our current paradigm for the formation of galactic struc- 
ture by gravitational instability in a theory with a primordial scale invariant 
spectrum of Gaussian perturbations. 

This work concentrates on the convergence field, which can be seen as a 
weighted projection of the density contrast field until the source. The character- 
istics of the lensing map will them be dependent not only on the density field, 
but on how it is weighted. I show how these two components relate to the cos- 
mic evolution, and what are the resulting lensing properties for two cosmological 
models. 

The differentiation of the lens-planes is small for the two models considered 
when the spectrum of matter fluctuations is normalized to match the abundance 
of rich clusters today. However the lensing geometrical factor, the term g/a, is 
clearly distinct between the two models, allowing the clear differentiation of 
the models by a power spectrum analysis, or by one-point statistics. It is to 
be expected, therefore, that for quintessence models these statistics will assume 
values intermediate between what was obtained for ACDM and OCDM, because 
the geometrical factor is intermediate. 

It would be very useful if the most popular N-body codes, for example Hydra 
(which was used here) , could incorporate quintessence models. That would allow 
a direct calculation of lensing measures in this model. 

The morphological analysis of the maps indicates some evolution of the 
second Minkowski functional (equal) for both models, and no distinguishable 
evolution of the third functional (genus) in the redshift range explored. The 
third functional is a better indicator of the departure from Gaussianity than 
the second functional. Because the third Minkowski functional depends weakly 
on the redshift of the source it is a specially recommended statistic to evaluate 
Gaussianity in surveys that have a poor redshift determination of the sources. 

References 

Barber A.J., 2002, astro-ph/0108273 
Bun E.F., White M. 1997, ApJ, 480, 6 

Carroll S.M., Press W.H., Turner E.L. 1992, ARA&A, 30, 499 
Couchman H.M.P., Thomas P. A., Pearce F.R. 1995, ApJ, 452, 797 
Guimaraes A.C.C. 2002, MNRAS, 337, 631 

Hamilton A.J.S., Kumar P., Lu E., Matthews A. 1991, ApJ, 374, LI 

Jain B., Scljak U., White S. 2000, ApJ, 530, 547 

Ma C.-P., Caldwell R.R., Bode P., Wang L. 1999, ApJ, 521, LI 

Scranton R., Dodelson S. 2000, preprint (astro-ph/0003034) 



